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Abstract 

The two-point correlation function of a Potts model on a graph 
G may be expressed in terms of the flow polynomials of 'Poissonian' 
random graphs derived from G by replacing each edge by a Poisson- 
distributed number of copies of itself. This fact extends to Potts 
models the so-called random-current expansion of the Ising model. 



1 Introduction 

The Tutte polynomial and its relatives have rarely been distant from the 
work of Dominic Welsh. They play important roles in matroid theory, [2U] . 
computational complexity, [221 I2H1 121] , and models of statistical physics, 
|211 l2lj . They provide the natural way to count and relate a variety of 
objects defined on graphs. We show here that they permit a representation 
of the two-point correlation function of a ferromagnetic Potts model on a 
graph G in terms of the flow polynomials of certain related random graphs. 
This representation extends to general Potts models the so-called random- 
current expansion for Ising models, wielded with great effect in [TJ |2l E3 UH] 
and elsewhere, and it amplifies the links between the Potts partition function 
and the Tutte polynomial surveyed earlier by Welsh and Merino, [21]. 

Two key elements of the analysis of the Ising model on a graph G are 
the random-cluster representation and the random-current expansion. The 
former is valid for all Potts models (and more besides), but the latter has not 
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previously been extended beyond the Ising model. It hinges on an expansion 
of the partition function in terms of 0/1-vectors indexed by edges and such 
that, for every vertex v, the sum of the values over edges incident to v is even. 
Such a vector may be recognised as a 'mod-2 flow'. It turns out that the q- 
state Potts partition function corresponds similarly to counts of 'mod-g flows' 
on a graph derived from G in the following way. Let A > 0, and replace every 
edge e of G by P(e) parallel edges, where the P(e) are independent Poisson- 
distributed random variables with parameter A. The quantity of interest is 
the mean number of non-zero mod-g flows on the resulting random graph. 

There is a powerful method of 'path-manipulation' by which many im- 
portant results have been proved for the Ising model. This method has a 
simple form when set in the context of a Poissonian random graph, and we 
illustrate this in Section 5 with a version of the 'switching lemma' of [1 . 

A short tour of graph polynomials appears in Section 2. In Section 3 
is introduced the Potts and random-cluster models, and the main result is 
proved in Section 4. Applications to the Ising model are summarized in 
Section 5. The principal open area is to extend the random-current analysis 
to Potts models with general q. 

2 Graph polynomials 

Let G = (V, E) be a finite graph, possibly containing multiple edges and 
loops. The Whitney and Tutte polynomials of G are well known to graph 
theorists, and we begin with a reminder of their definitions. The (Whitney) 
rank- generating function of G was introduced in [25^ and is given by 

W G (u,v) = ^ u r(G 'V (G '>, u,veR, (2.1) 

E'CE 

where r(G') = \V\ — k(G') is the rank of the subgraph G' = (V,E r ), and 
c(G') = \E'\ — \V\ + k(G') is its co-rank. Here, k(G') denotes the number of 
components of G' . Note that 

W G (u,v) = (u/v) lvl v ]E ' l (v/u) k(E '\ u 7 v^0. (2.2) 

E'CE 

The rank-generating function has various useful properties, and it occurs 
in several contexts in graph theory, see [HIES]- The Tutte (or dichromatic) 
polynomial of G was introduced independently in ^HJ El] , and may be ex- 
pressed as 

T G (u, v) = (u- l^-'Wcdu - v-1). (2.3) 
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This also is a function of two variables. For suitable values of these variables, 
it provides counts of colourings, forests, and flows, and of other combinatorial 
quantities. The principal purpose of the current paper is to explore the use 
of the Whitney/Tutte polynomial in the study of the correlation functions of 
the Potts model, and to this end we define next the flow polynomial of G. 

We turn G into a oriented graph by allocating an orientation to each 
edge e £ E, and we denote the resulting digraph by G = (V,E). If the 
edge e = (u, v) £ E is oriented from u to v, we say that / leaves u and 
arrives at v. It will turn out that the choices of orientations are immaterial 
to the principal conclusions that follow. Let q £ {2, 3, . . . }. A function 
/ : E — > {0, 1, 2, . . . , q — 1} is called a mod-q flow on G if 

2^ f{z) ~ /J /(e) = modulo q, for all v £ V, 

eel: e £ E : 

e leaves v § arrives at v 

which is to say that flow is conserved at every vertex. A mod-g flow / is 
called non-zero if /(e) ^ for all e £ E. Let Cc{q) be the number of non- 
zero mod-g flows on G. It is fundamental that the quantity Cg{q) does not 
depend on the orientations of the edges of G, and the proof may be found 
in [H]. The function Cc{q), viewed as a function of q, is called the flow 
polynomial of G. 

The flow polynomial may be obtained as an evaluation of the Whit- 
ney/Tutte polynomial with two particular parameter values, as follows: 

C G (q) = (-l)^W G (-l,-q) 

= (-1)1*1-1^1+1^(0, 1 _ q ) 7 q G {2 , 3, . . . }. (2.4) 

See EH]- We shall later write C{G\q) for Cc(q), and similarly for other 
polynomials when the notational need arises. 



3 Potts and random-cluster models 

Amongst models for ferromagnetism, the Potts model is one of the most 
studied. It has two principal parameters, the 'inverse temperature' (3 £ 
(0, oo) and the number q £ {2, 3, . . . } of local states. When q = 2, the Potts 
model becomes the Ising model. Let G = (V, E) be a finite graph which for 
simplicity we assume to have no loops. It is convenient to allow a separate 
parameter for each edge of G, and thus we let J = (J e : e £ E) be a vector 
of non-negative numbers, and we set 

p e = 1 - e _/3Je9 , eeE. (3.1) 



3 



The configuration space of the g-state Potts model on G is the set S 
{1,2,..., q} v . The Potts measure on £ is given by 

7T/3J,g(o-) = -^p expl^ 0J e (q5 e (a) - 1)1, a G E, (3.2) 
where, for e = (x, y) G i£, 



5 e (a) = b ax ^ y 



1 if o"z = <? yi 
otherwise, 



and Z p = is the partition function 

Z p = ^exp|^/3J e (g5 e (^) - 1)}. (3.3) 

Since (3J e > 0, the Potts measure 7Tpj )q allocates greater probability to con- 
figurations for which 5 e (a) = 1 for a larger set of edges e. That is, it prefers 
configurations in which many neighbour-pairs have the same state, and in 
this regard the model is termed 'ferromagnetic'. 

A central quantity is the 'two-point correlation function' given by 

TpJ*(x, V) = TT/3J,g{^x = Vy) ~ ~, X, y G V. (3.4) 

We shall work here with gTgj i(? (x, y) = ^i3j,q(qS ax ^ y — 1) and, for ease of 
notation in the following, we write 

a(x, y) = qrp J>q (x, y), x, y G V, (3.5) 

thereby suppressing reference to the parameters (33 and q. 

Two of the most successful ways of studying the Ising/Potts models are 
the so-called 'random-cluster model' and the 'random- current expansion'. 
We define next the random-cluster model, and we explain its relevance to 
the Potts model. The random-current expansion for the Ising model will be 
reviewed in Sectional 

In the (bond) percolation model on G, each edge is declared at random 
to be either 'open' or 'closed'. An edge is declared 'open' with some given 
probability p, and closed otherwise, and different edges are allocated inde- 
pendent states. The percolation model is basic to the study of disordered 
media, particularly when the underlying graph is part of a 'crystalline' lattice 
such as the d-dimensional cubic lattice h d . See E3 for a full account. When G 
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is a complete graph, the percolation model is usually called an 'Erdos-Renyi 
random graph', see [T3j . 

The random-cluster measure on G is obtained through a perturbation 
of the percolation measure, as follows. Let p = (p e : e G E) G [0, 1} E and 
q G (0, oo). The configuration space is Q = {0, 1} E . For ueO and e 6 £, we 
say that e is u-open (or, simply, open) if u;(e) = 1, and u-closed otherwise. 
The random-cluster probability measure on Q is defined by 




(i)Gfl, 



where k(u) denotes the number of cu-open components on the vertex-set V, 
and Z RC = Zq C is the appropriate normalizing factor. We sometimes write 
4>G,p,q when the role of G is to be emphasized. 

It is common to take p e = p for all e G E, in which case we write for 
0p ; g. The special case q — 1, p = p is evidently the percolation measure with 
parameter p, in which case we write <j) p = <p Pi i. It turns out that the random- 
cluster model with q G {2, 3, . . . } corresponds in a certain way to the Potts 
model on G with q local states and with ft J satisfying (j3.1|) . Specifically, the 
two-point correlation function of the latter is (up to a harmless factor) equal 
to the connection probability of the former, 

Tpj, q (x, y) = (1 - q- l )<p^ q (x «-> y), x, y G V, (3.6) 

where x <-> y means that there exists a path of open edges from x to y. 
The random-cluster model was introduced by Fortuin and Kasteleyn around 
1970, and has been reviewed recently in jHHHHHI]. 

The random-cluster partition function Zq C is given by 

Z% c (p,q) = J2p M " )l (l-p) lE \ vHl q HuJ \ 
and is easily seen by (|2.2|) to satisfy 

z? M = - P rw c (-^— y , P/ 1, 

a relationship which provides a link with other classical graph-theoretic quan- 
tities. See 00 0III1 12U- 

4 Potts correlations and flow counts 

It is shown in this section that the Potts correlation functions ()3.4)1 may be 
expressed in terms of flow polynomials associated with a certain 'Poissonian' 
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random graph derived from G by replacing each edge by a random number of 
copies. This extends to general q the random-current expansion of the Ising 
model described in Sectional 

For any vector m = (m(e) : e G E) of non-negative integers, let G m = 
(V, E m ) be the graph with vertex set V and, for each e G E, with exactly 
m(e) edges in parallel joining the endvertices of the edge e; the original edge 
e is itself removed. Note that 



|£ m | = ^m(e). (4.1) 



eG-B 

Let A = (A e : e G E) G [0, oo) E . Let P = (P(e) : e G £) be a family of 
independent random variables such that P(e) has the Poisson distribution 
with parameter A e . The random graph Gp = (V,Ep) is called a Poisson 
graph with intensity A. Let and denote the corresponding probability 
measure and expectation operator. 

For x, y G V, x ^ y, we denote by G F ,y the graph obtained from Gp by 
adding an edge with endvertices x, y. If x and y are already adjacent in Gj>, 
we add exactly one further edge between them. Potts-correlations and flows 
are related by the following theorem. The function a(x,y) is given in ()3.5|) . 

Theorem 4.2. Let q G {2, 3, . . . } and X e = (3J e . Then 

E x (C(G^;q)) 

a ^ y)= E x (C(G P ;q)) > X > yEV - (4 - 3) 

This formula is particularly striking when q = 2, since non-zero mod-2 
flows necessarily take only the value 1. A finite graph H = (W, F) is called 
even if the degree of every vertex w is even. It is trivial that C^(2) = 1 if if 
is even, and C#(2) = otherwise. By ()4.3j) . for any graph G, 

.fev) = P ; ( ff' g . iSeVen , ) . (4.4) 
Pa(G p is even) 

Such observations are at the heart of the random- current expansion for Ising 
models. See Sectional 

Theorem 14.21 may be extended via ()3.6J) to the random-cluster model. 
Assume for simplicity that every edge has the same parameter p. The proof 
of the following is easily derived from Theorem 14 .2^ and may be found in 
[llj. It is obtained by expressing the flow polynomial in terms of the Tutte 
polynomial T, and allowing q to vary continuously. 

Theorem 4.5. Let p G [0,1) and q G (0, oo). Let A e = A for all e G E, 

where p = 1 — e~ Xq . 
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(i) For x, y G V , 



( 1U ( \ E A ((-l)^lT(^;0,l-g)) 

(g - - y) = EA(( _ 1)W(Gp;0>1 _ g)) • ( 4 - 6 ) 

(ii) Forge {2,3,...}, 

p (g*< w) ) = (1 -p)^- 2 ^q^E x (C(G F ;q)). (4.7) 

When q = 2, ()4.7|1 reduces to the curiosity 

p (2 fc(w) ) = 2 |y| P A (G P is even). (4.8) 

This may be simplified further. Let ((e) = P(e) modulo 2. It is easily seen 
that Gp is an even graph if and only if is even, and that the C( e ); e G E, 
are independent Bernoulli variables with 

PA(C(e) = l) = !(l-e- 2A ) = b. 
Equation ()4.7j) may therefore be written as 

(t) p (2 k{u) ) = 2 |v V p/2 (the open graph on V is even). (4.9) 



Proof of Theorem \4-£\ Since the parameter (3 appears together with the mul- 
tiplicative factor J e , we may without loss of generality take (3=1. We begin 
with a calculation involving the Potts partition function Z p of ()3.3|) . Let 
Z + = {0, 1, 2, ... } and let m = (m e : e G E) G Z+. By a Taylor expansion 
in the variables J P , 



exp 



Jin, 



J=0 



e€E > m( zzf V eG£ 7 

E A (9 P Z P | ) (4.10) 



where 
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By flU with p = \, 



d m z P 



3=0 



= e n - x ) 

crgE e€£ m 

= E [] £ h 5 n ei + 5n e ,l^e(0-)] 
crgS egS m n e S{0,l} 

= ^ ^(_l)|{e:»e=0}] ? |{e:ne=l}| f JJ 
ne{0,l} E ™ o-GE \&Em 

— ^_^y{ e; ne=0}lgl{ e: ™e=l}|^fc(m,n)^ (4-H) 

ne{0,l} B "i 

where k(m, n) is the number of connected components of the graph obtained 
from G m after deletion of all edges e with n e = 0. By (J2.2j) - (j2.4|) . 



d m Z p 



— (_l)l S m| fq\\W-ne=l}\qk(m,n) 

j=o ' 

ne{0,l} B "i 



= (-l)\^\ q W\W Gm (-l,-q) (4.12) 

= g |V ^(G m ;g). (4.13) 

Combining (OTlll - lHl^l . 

expj- ^ J e U p = gl y lE A (C(G P ; g)). (4.14) 

^ e£E ' 

Let x,y (z V. We define the unordered pair / = (x, y), and write Sf(a) = 
5 axtay for cr G S. We have that 



cr(x,y) = npj jq (q5 f (a) - 1) 

= Jp E^ ff ) " l)exp{X;^(?W - 1)}. 



(4.15) 



By an analysis parallel to (|4~TD|l - (|4~T4]l . 

exp[-^ A J2(q6 f {<r) - l)exJj2P-UqW ~ 1)} 

^ e£E > cr£E ^e£E ' 

= q\ y ^ x {C{G^-q)), 
and (jOJ) follows by (ETHl and (05]) . □ 



(4.16) 
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5 Random- current expansion of the Ising model 



Unlike the situation with the Potts model, there is a fairly complete analysis 
of the Ising model. A principal part in this analysis is played by Theorem 
14.21 with q = 2, under the heading 'random-current expansion'. This has per- 
mitted proofs amongst other things of the exponential decay of correlations 
in the low-/? regime on the cubic lattice L d with d > 2. See 121 E|. It has 
not so far been possible to extend this work to general Potts models, but 
Theorem 14.21 could play a part in such an extension. 

Let G = (V, E) be a finite graph without loops as before, and set q = 2. 
We restrict ourselves here to the Ising model with J e = J for all e G E, and 
we write A = (3 J. By Theorem 14. 2\ 

a(x, y) = 2r A , 2 (x, y) = ^ even) o < A < oo. (5.1) 

Pa(Gp is even) 

The value of such a representation will become clear during the following 
discussion, which is based on material in QjJ EE] ■ in advance of this, we 
make a remark concerning (j5.1|) . In deciding whether Gp or Gp V is an even 
graph, we need only know the numbers P(e) when reduced modulo 2. That 
is, we can work with ( G Q = {0, 1} E given by ((e) = P(e) mod 2. Since 
P(e) has the Poisson distribution with parameter A, ((e) has the Bernoulli 
distribution with parameter 

p> = P A (P( e ) is odd) = f (1 - e" 2A ). 

We obtain thus from (|5.1|) that 

_ <P p ,(d( = {x,y}) 

{ ,y) ~ Md( = z) ' 

where </y denotes product measure on Q with density p', and 
d( = LeV: C(e) is oddj, (en, 

^ e: e^v ^ 

where the sum is over all edges e incident to v. We refer to members of d( 
as 'sources' of the configuration (. 

Let M = (M e : e G P) be a sequence of disjoint finite sets (possibly 
empty) indexed by E, and let m e = \M e \. As noted in the last section, the 
vector M may be used to construct a multigraph G m = (V, E m ) in which 
each e G E is replaced by m e edges in parallel; we may take M e to be the set 
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of such edges. For x, y G V, we write l x <-» y in m' if x and y lie in the same 
component of G m . We define the set <9M of sources of M by 

8M=\veV: ^ m e is odd I. (5.2) 

Thus, for example, G m is even if and only if <9M is empty. From the vector 
M we construct a vector N = (N e : e G E) by deleting each member of each 
M e with probability |, independently of all other elements. That is, we let 
Bi-, i G lJ e M e , be independent Bernoulli random variables with parameter 
~, and we set 

N e = {iE M e : Bi = 1}, e G £. 

We write P M for the appropriate probability measure. The following lemma 
is pivotal for the computations which follow. 

Lemma 5.3. Let M and m be as above. Ifx, y G V are such that x ^ y and 
x «-> y in m £/ien ; /or A C V , 

P M (9N = <9(M\N) = A) = P M (<9N = 0, d(M \ N) = A A 

Proof. Take M e to be the set of edges of G m parallel to e, and assume that 
x <-> ?/ in m. Fix A C V. Let be the set of all vectors n = (n e : e 6 £) 
with n e C M e for all e. Let p be a fixed path of G m with endpoints x, y, 
and consider the map p : .M — >• M. given by 

p(n) = n A p, n G AL 

The map p is one-one, and maps {n G .M : <9n = 9(M \ n) = A} 

to {n' G : 9n' = 0, <9(M \ n') = A A {a;,y}}. Each member of M. is 
equiprobable under P M , and the claim follows. □ 

Let A G [0, oo), and recall from the last section the definition of a Poisson 
graph with parameter A. The following is a fairly immediate corollary of the 
last theorem. Let M = (M e : e G E) and M' = (M' e : e G E) be vectors 
of disjoint finite sets satisfying M e n Mj. = for all e, / G E, and suppose 
that the random variables m e = \M e \, m' e = \M' e \, e G E are independent 
and such that, for each e G E, m e and mj, have the Poisson distribution 
with parameter A. We write MUM'= (M e U M' e : e G £), and P for the 
appropriate probability measure. The following lemma is a simplification of 
the so-called switching lemma of pQ. 
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Lemma 5.4. If x,y G V are such that x ^ y and x *->■ y inm + m' then, 
for ACV, 

P(<9M = {x, y}, dM' = A | MUM') 

= F(dM = 0, dM' = A A {x, y} | M U M') P-a.s. (5.5) 

Proof. Conditional on the sets M e U M' e for e G E, the sets M e are selected 
by the independent removal of each element with probability ~. The claim 
now follows from Lemma f5. HI □ 

We present two applications of Lemma 15.41 to the Ising model, as in . 
For m = (m e : e G E) G Zf , let 



9m = |t> G : m e is oddj-. 



(5.6) 



as in (|5.2j) . We write as before 

a(x, y) = 2r Ai2 (x, y) = 7r Aj2 (25 CT:ciCTy - 1), x, y G V, 

thereby suppressing reference to A. By (15. lj) . 

P A Q9P = {x,y}) 

cr(x,v) = — — 77— : — . (5.7) 

V ' y; P A (dP = 0) V 7 

Let Qa denote the law of P conditional on the event {<9P = A}, 

Q A (E) = F x (PeE\dP = A). 

We have need of two independent copies Pi, P2 of P with potentially different 
conditionings, and thus we write Qa-b = Qa x Qb- 

Lemma 5.8. Let x,y,z G V be distinct vertices. Then: 

(i) <t(x, yf = Q 0]0 ( x <-> V in p i + p 2), 

(ii) a(x, y)a(y, z) = a(x, z)Q{ XjZ y. (x <-> y in P 1 + P 2 ) . 

Proof, (i) By ()5.7|) and Theorem 15. 4[ 



2 P A x P A (9Px = {x, y}, <9P 2 = {ar, y}) 



P A (dP = 0) 2 
P A x P A (gPi = {x, y}, dP 2 = {x, y},x^y in Pj 

P A (<9P = 0) 2 
P A x P A (9Pi = <9P 2 = 0, x <-> y in P x + P 2 ) 

P A (<9P = 0) 2 
Q 0;0 (x^y in P1 + P2). 
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(ii) Similarly, 

<r(x,y)<r(y,z) 

_ Pa x P A (dPi = {x,y}, dP 2 = {y,z}) 

P A (<9P = 0) 2 

_ P A x P A (<9Pi = 0, <9P 2 = {x, z}, x «-> y in P x + P 2 ) 
= P A (<9P = 0) 2 

= PA lf P ' = {g ' * }) • Pa x P A (x <-> y in P, + P 2 | 9P X = 0, dP 2 = {x,z}) 
\r\(oP = 0) v 

= cr(x, z)Q{ x ,z};0(a; <-> y in Pi + P 2 ), 

and the proof is complete. □ 

Theorem (|5.8j) (ii) leads to an important correlation inequality known as 
the 'Simon inequality', [IB]. Let x,z G V be distinct vertices. A subset 
W C V is said to separate x and z if x, z W and every path from x to z 
contains some vertex of W. 

Theorem 5.9. Let x,z £ V be distinct vertices, and let W separate x and 
z. Then 

a(x,z) < ^(x,y)(r(y,z). 
yew 

Proof. By Theorem 15.8( 11). 

o-(x,y)a(y,z 



— — = 2^ ®{**yA x H f mp i + p 2) 

V 

x,z}-,0{\{y e W : x «-> y in P : + P 2 }|). 



Assume that the event dPi = {x, z} occurs. On this event, x «-> 2 in Pi + P 2 . 
Since separates x and z, the set {y G : x «-> y in P! + P 2 } is non-empty 
on this event. Thus its (conditional) mean size is at least one under Q{ Xt z};0, 
and the claim follows. □ 

The Ising model on the graph G = (V, E) corresponds as described in 
Section El to a random-cluster measure 4>p,q with q = 2. By (|3.6|) . 

cr(x, y) = 2r A , 2 (x, y) = (f) Pi g(x y), 

where p = 1 — e _Al3 and g = 2. The Simon inequality may be written in the 
form 

0p,g(x <-> 2) < Y &P>v( X *~* V^liV *-* Z ) 
ydW 
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whenever W separates x and z. It is a curious fact that this inequality holds 
also when q — 1, see j§] Q3- One ma Y conjecture that it holds for any 
qe[l,2\. 

Let d > 2. The random-cluster measure Pi9 on the cubic lattice h d 
may be obtained as a weak limit (with so-called free boundary conditions) 
of the random-cluster measure on finite boxes A, as A f Z d . The percolation 
probability is the function 9 given by 

9(p,q) = (j) m (0 *-> oo), 

the probability that the origin is the endpoint of an infinite open path. The 
critical point is defined as 

p c (q) = sup{p : 9(p, q) = 0}. 

Let || ■ || be a norm on Z d . It has been conjectured that, for p < p c {q), 
there exists 7 = j(p,q) G (0, 00) such that 4> Pt q(0 <-> 2) < e - "^" 7 for all 
x E 1< d . This has been proved when q — 1,2 and g is sufficiently large. The 
Simon inequality implies the following necessary and sufficient condition for 
exponential decay when q = 2. 

Theorem 5.10. Let q = 2 and assume that p is such that 

0p,g(O <-> X) < OO. 

T/iere exists 7 = 7(p, g) G (0, 00) s«c/i 

<f>p, q {0 ^ z)< e" l|z|17 , 2 G Z d . 

The proof follows standard lines, and may be found in QBj together 
with proofs of the following facts. There is an important extension of the 
Simon inequality due to Lieb, [IS]. This also may be proved via the flow 
representation of Theorem 14.21 The Lieb inequality has an important conse- 
quence for the nature of the phase transition of the Ising model, namely the 
'vanishing of the mass gap'. 

Let q = 2 and write 

V>(p, q) = lim I -- log Pi9 (O <-> <9A„) 1 , 

n^oo [ f| J 

where A n = [— n, n] d and dA n = A n \A n _i. Note that ijj(p, q) is non-increasing 
in p, and 4>(p,q) = if p > p c (q)- One of the characteristics of a first-order 
phase transition is the (strict) exponential decay of connectivity probabilities 
at the critical point. The quantity ipipdq)^) is sometimes termed the mass 
gap. 
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Theorem 5.11. Let d > 2 and q = 2. Then ip(p,q) decreases to as 
P T Pc(q)- I n particular ip(p c (q), q) = 0, that is, the mass gap equals 0. 

This was proved in see also [TT]. The corresponding statement is 
known to be false for d > 2 and q > Q(d) for some sufficiently large Q(d). 
See jTHIHj. In further use of the random- current expansion (with q = 2), it 
has been proved that i/;(p, q) > whenever p < p c (q). See 13 El for more 
details of the Ising phase transition. 
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